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ABSTRACT 

i-£h ■ Context. Studying molecular gas in the central regions of the star burst galaxies NGC 4945 and Circinus enables us to 

characterize the physical conditions and compare them to previous local and high-z studies. 

Aims. We estimate temperature, molecular density and column densities of CO and atomic carbon. Using model pre- 
» I . dictions we give a range of estimated CO/C abundance ratios. 

Methods. Using the new NANTEN2 4m sub-millimeter telescope in Pampa La Bola, Chile, we observed for the first 
time CO 4-3 and [C i] 3 Pi — 3 Po at the centers of both galaxies at linear scale of 682 pc and 732 pc respectively. We 
compute the cooling curves of 12 CO and 13 CO using radiative transfer models and estimate the physical conditions of 
CO and [CI]. 

Results. The centers of NGC 4945 and Circinus are very [C i] bright objects, exhibiting [C i] 3 Pi - 3 Po luminosities 
of 91 and 67Kkms _1 kpc 2 , respectively. The [Cl] 3 Pi — Po/CO 4-3 ratio of integrated intensities are large at 1.2 in 
! NGC 4945 and 2.8 in Circinus. Combining previous CO J= 1-0 , 2-1 and 3-2 and 13 CO J= 1-0 , 2-1 studies with 

. our new observations, the radiative transfer calculations give a range of densities, n(H2) = 10 3 — 3 x 10 4 cm -3 , and a 

0^ ' wide range of kinetic temperatures, Tki n = 20 — 100K, depending on the density. To discuss the degeneracy in density 

and temperature, we study two representative solutions. In both galaxies the estimated total [CI] cooling intensity is 
stronger by factors of ~ 1 — 3 compared to the total CO cooling intensity. The CO/C abundance ratios are 0.2-2, similar 
to values found in Galactic translucent clouds. 

Conclusions. Our new observations enable us to further constrain the excitation conditions and estimate the line emis- 
sion of higher- J CO- and the upper [CI]-lines. For the first time we give estimates for the CO/C abundance ratio 
in the center regions of these galaxies. Future CO J= 7-6 and [CI] 2-1 observations will be important to resolve the 
ambiguity in the physical conditions and confirm the model predictions. 



Key words. Keywords should be given 

1. Introduction Emission of CO, [Cl] (and [Cn]) traces the bulk of car- 
bon bearing species in molecular clouds that play an impor- 
The spiral galaxies NGC 4945 and Circinus at distances tant role in their chemical network. The CO 4-3 transition 
of ~ 3.7 and ~ 4Mpc belong to the nearest and infrared in particular is a sensitive diagnostic of the dense and warm 
brightest galaxies in the sky. Their strong central star burst gas while the CO 1-0 transition traces the total molecular 
activity is fed by large amounts of molecular material and mass. The [Cl] 3 Pi — 3 Pq (henceforth 1-0) line has, to date, 
this has been studied extensively at millime ter wavelengths been detected in about 30 galactic nu clei (jGerin fc Phillipsl 
(e.g. lCurran et al.ll2001aHWang et al.ll2004f) . However, sub- [2000l ; llsrael feB aas 200^ llsraeil |2005|) and appears to trace 
millimeter observations are largely missing. The important the surface regions of clumps ill uminated by the FUV fiel d 
rotational transitions of CO with J > 4 and the fine struc- of newborn, massive stars (e.g. iKramer et al.ll2004l . |2005[ ). 
ture transitions of atomic carbon have not yet been ob- The use of CI as an accurate trace r of the cloud mass 
served. These transitions often contribute significantly to has been discussed with controversy dFrerking et al.lll989l ; 
the thermal budget of the interstellar gas in galactic nuclei iPapadopoulos fc Grevell2004t iMookeriea et al.ll2006f ). 
and are therefore important tracers of the physical condi- 
tions of the warm and dense gas. The strong cooling emission is balanced by equally 

strong heating caused by the vigorous star formation ac- 

tivity in the galaxy centers. The variation of CO cooling 

Send offprint requests to: M. Hitschfeld, e-mail: intensities with rotational number, i.e. the pea k of the CO 

hitschfeld@phl.uni-koeln.de cooling curve, reflects the star forming activity (jBavet et al.l 
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I2006T) and, possibly, also the underlying heating mecha- 
nisms. Several mechanisms have been proposed to explain 
the heating of the ISM in galactic nuclei and it is cur- 
rently rathe r unclear which of these dominates in individual 
sources (e.g. lWang et~aHl2004l) . Heating by X -rays from the 
active galactic nuclei (AGN) may lea d to strongly enhance d 
intensities of high- J CO transitions (M ciierink et alJ[2007t ). 
possibly allowing one to discriminate this heating mecha- 
nism from e.g. stellar ultra violet heating via photon domi- 
nated regions (PDRs) (e.g. lBavet et al.ll2006[ ). The greatly 
enhanced supernova rate by several orders of magnitude 
relative to the solar system value leads to an enhanced 
cosmic ray flu x, providing another source of gas heatin g 
in the centers (|Farquhar et al.lll994t [Bradford et al.ll2003f) . 
Another mechanism for heating is provided by shocks. The 
most important, large-scale shocks are produced by den- 
sity wave instabilities which induce gravitational torques 
(in spiral arms and/or bars ) and make the gas fall into the 
nucleus (jUsero et al.ll2006f ). A non-negligible contribution 
is also given by shocks produced by supcrnovac explosions. 
On smaller scales, bipolar outflows from young stellar ob- 
jecs (YSO) can also contribute to this heating , although to 
a lesser extent (e.g. iGarcia- Burillo et al.ll2001l) . 



obtained for NGC4945. Large amounts of molecular gas 
have been found by s tudies of low-J CO observations 
(I Johansson et all 1991 jAalto et all 119951 : lElmouttie et al.l 
119981: iCurran et al J 1 1998ft including C 17 Q. C 18 and HCN 



(ICurran et al.l l2001af). The stron g H2O maser emission 



found ( Gardner fe Whiteoaklll982l ) traces a thin accretion 
disk of 0.8 pc radius, with a significant population of 
masers lying away fro m this disk, possibly in an outflow 
(jGreenhill et al.lll99 8). Its obscured nucleus is classified as 
an X-ray Compton thick Seyfert 2. It sho ws a circumnu- 
clear star burst on scales of 100 - 200 pc (jMaiolino et al.l 
Il998ft with a complex structure of H 11 as seen in Hubbl e 
Space Telescope (HST) Ha images (| Wilson et al.ll2000D . 
Adaptive optics studies find there has been a recent star 
burst (~ 100 Myr old) in t he central 8pc accounting fo r 
2% of the total luminosity (|Mueller Sanchez et afll2006ft . 
The strong FIR emission, the large molecular gas reser- 
voir, the existe nce of a molecular ring associated with star 
burst activity (jCurran et al.lfl998l ) and the similarity with 
NGC 4945 make them ideal objects for comparative studies 
of the dense and warm ISM in their nuclei. 

2. Observations 



1.1. NGC 4945 

NGC 4945, a member of the Centaurus group of galaxies, 
is seen n early edge-on (Table[I| with an optical diameter 
of ~ 20' (|de Vaucouleurs et al.lll99i[ ). Hi kinematics indi- 



cate a galaxy mass of 1.410 11 Mq within a radius of 6.3' 
with molecul ar and neutral atomic gas contributing ~ 2% 
respectively (|Ott et al.ll200~l[) . 

Wit h a dynamical mass of ~ 3 10 9 M Q in the central 
600 pc (jMauersberger et aT1ll996ft , it is one of the strongest 
IRAS point sources with almost all the far infrared lu - 
minosity coming from the nucleus (|Brock et alj Il988l ). 
Observations of the X-ray spectrum are consistent with a 
Seyfert nucleus (jlwasawa et al.|[l993[ ). and further analysis 
of op tical imaging and infrared spectra (jMoorwood fe Olival 
|1994|) suggests that this object is in a late stage on the tran- 
sition from starburst to a Seyfert galaxy. Its nucleus was 
the first s ource in which a powerful H2O mega maser was 
detected (|Dos Santos fc Lepindll979ft . 



Stud ies in Hi (lAbles et all 1 1987ft and low-J CO t ran- 
sitions dWhitcoa k et all 11990: IDahlem et al.l Il993t " 



119951 : IMauersberger et ahl 11996ft suggest the presence 



of a face-on circumnuclear molecula r ring. Mill i meter 
molecular multiple transition studies dWang et al.l 120041 : 
ICunningham fc Whiteoakll2005ft are consistent with this re- 
sult. The bright in frared and radio emission in the nucleus 
(jGhosh et al.lll992ft , and the evidence t hat large amounts o f 
gas seem to coexist in the central 30" (|Henkel et al.lll994ft 
make it particularly suited for studying the high density 
environment in the center of this galaxy. 



1.2. Circinus 

The nearby starburst spiral Circinus has a small optical 
angular diameter of ~ 7'compared to its hydrogen diame- 
ter £> H = 36' defined by the 1 10 20 atoms cm 2 contour in 
iFreeman etal] (|1977l ). 

Circinus has a dynamical mass of ~ 3 10 9 M Q within 
the inner 560 pc (|Curran et al.l [i"998ft matching the value 



Table 1 . Basic properties of Circinus, NGC 4945 and NGC 253. 
LiR is the total mid and far infrared luminosity calculated 
from the flux densities at 12/im, 25/^m, 60/zm, 100pm listed in 
the IRAS point source catalogue (|Lonsdale fc H clou 1985]) us- 
ing the formulae given in Table 1 of iSanders fc Mirabell (119961 ) 
and the distances listed below. References: a RC3 catalogue 
o f Ide Vaucouleurs et all (| 19911), b IMauersberger et al.l ([19961), 
c IFreeman etHI (fl977l ) d ICurran et~aTl (jl998l), e iRekola et all 
(120051 ), 1 ICurran et all (|2001bl ). g lFullmer & Londsdald (|1989l ) . 





Circinus 


NGC 4945 


NGC 253 


RA(2000) 


14:13:09.9 


13:05:27.4 


00:47:33.1 


DEC(2000) 


-65:20:21 


-49:28:05 


-25:17:18 


Type 


SA(s)b a 


SB(s)cd a 


SAB(s)c a 


Distance [Mpc] 


4.0 d 


3.7 6 


3.5 e 


38" correspond to 


732 pc 


682 pc 


646 pc 


LSR velocity [kms~ ] 


434 


555 


243 


Inclination [deg] 


65 c 


78 6 


78 a 


Lm [10 10 L ] 


1.41 9 


1.39 9 


2.67 9 


5ioo[Jy] 


3.16 10 2 


6.86 10 2 


10.4 10 2 



We used the new N ANTEN2 sub-millimeter telescope (e.g. 
iKramer et aT1l2007ft on Pampa la Bola at an elevation of 
4900m together with a dual channel 490/810 GHz receiver 
(operated jointly by the Nagoya University radioastron- 
omy group and the KOSMA group from Universitat zu 
Koln) to observe the centers of NGC 4945, Circinus, and 
NGC 253 (Table [J) in CO 4-3 and [Ci] 1-0. The obser- 
vations were performed from September to October 2006 
using the position-switch mode with 20 sec On- and 20 sec 
Off-time. In October 2007 we reobserved 12 CO 4-3 in the 
centers of NGC4945 and Circinus with the newly avail- 
able chopping tertiary in double beam-switch mode yielding 
significantly improved baselines and removing atmospheric 
features in the spectra. In NGC 4945 the total integration 
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times ON-source are 17min and lOmin for 12 CO 4-3 and 
[C i] 1-0 respectively and in Circinus integration times ON- 
source are lOmin and 34min for 12 CO 4-3 and [Ci] 1-0 
respectively at source elevations of 40-60°. The relative cal- 
ibration uncertainty derived from repeated pointings on the 
nuclei is about 15%. 

Position-switching was conducted by moving the tele- 
scope 10' in azimuth, i.e. out of the galaxy. Using dou- 
ble beam-switch (dbs) mode, the chopper throw is fixed 
at 162" in azimuth with a chopping frequency of 1 Hz. 
Typical double-sideband receiver temperatures of the dual- 
channel 460/810 GHz receiver were ~250K at 460 GHz and 
492 GHz. The system temperature varied between ~850 
and ~1200K. As backends we used two acusto optical spec- 
trometers (AOS) with a bandwidth of 1 GHz and a chan- 
nel resolution of 0.37 kms -1 at 460 GHz and 0.21 kms -1 at 
806 GHz. The pointing was regularly checked and pointing 
accuracy was stable with corrections of ~ 10". The half 
power beam width (HPBW) at 460 GHz and 492 GHz is 
38" with a beam e fficiency B e g = 0. 50 and a forward effi- 
ciency F ff = 0.86 (|Simon et al.ll2007t ). The calibrated data 
on T^-scale were converted to T m b by multiplying by the 
ratio F e g/B e g. The standard calibration procedure derives 
the atmospheric transmission (averaged over the bandpass) 
from the observed difference spectrum of hot load and blank 
sky, the latter taken at a reference position. Next, the model 
atmosphere atm is used to derive the atmospheric opacity 
taking into account the sideband imbalance. This is an im- 
portant correction, especially when observing the CO 4-3 
and [CI] 1-0 lines. This pipeline produces the standard spec- 
tra on the antenna temperature scale (TJ). In addition, we 
removed baselines up to first order. 

The 810 GHz channel was not used for these observa- 
tions due to insufficient baseline stability. 

All data presented in this paper are on the T m b scale. 

2.1. 12 CO 4-3 in NGC253 

To check the NANTEN2 calibration scheme, we re- 
trieved the 12 CO 4 -3 map of NGC253 taken at APEX 
(jGiisten et al.ll2006[ ) and smoothed it to the NANTEN2 
HPBW of 38" using a Gaussian kernel. The resulting spec- 
tra at the center position are shown in Fig.[T] Line temper- 
atures and shapes show very good agreement. 
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Velocity (km/ s) 

Fig. 1. Observed 12 CO 4-3 line emission in NGC253 from 
APEX (black line) and NANTEN2 (red line). We attribute 
small but systematic deviations of peak temperatures in the 
blue wing of the 12 CO 4-3 line to small pointing variations. 
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3. Spectra of CO 4-3 and [C i] 1-0 

Figured shows the 12 CO 4-3 and [Ci] 1-0 spectra of 
NGC 4945 and Circinus obtained with the NANTEN2 tele- 
scope. CO 4-3 spectra peak at 700 mK in NGC 4945 and 
250 mK in Circinus. Outside the velocity ranges of 350- 
800 kms -1 and 200-600 kms -1 respectively, the baseline 
rms values are llmK and 25 mK, respectively, at the ve- 
locity resolution of 15 kms" 1 . [Ci] 1-0 spectra peak at 
900 mK in both galaxies while the rms values are 110 mK 
and 140 mK, respectively. The [Ci] 1-0 area-integrated 
luminosities are 91Kkms _1 kpc 2 and 67Kkms _1 kpc 2 in 
N GC4945 and C ir cinus. 

ICurran et all (|2001al 1 19981) and iMauersberger et all 
(|l996l ) mapped the low- J 2-1, and 3-2 transitions of CO 
in the centers of both galaxies with SEST. We smoothed 
these data to the resolution of the NANTEN2 data, i.e. to 
38". Only the 12 CO 3-2 spectrum in Circinus is shown at 
its original resolution of 15" because a map of the central 
region could not be retrieved. Calibration of the CO 3-2 
spectrum in NGC 49 45 was confirmed recently at APEX 
dRisacher et al.ll2006l) . The 12 CO 1-0, 13 CO 1-0 and 13 CO 
2-1 spectra of th e central region of NG C 4945 and Circinus 
are presented in ICurran et al.l (|2001a| ). We list integrated 
intensities in Tabled] 

NGC 4945 shows broad emission between 350 and 
800 kms" 1 . The CO 4-3 line shape resembles the line shapes 
of 1-0 and 2-1. The line shape of [C i] 1-0 is similar to that 
of the CO transitions with a slighlty higher peak tempera- 
ture than CO 4-3. 

In CO 1-0 and 2-1, Circinus shows broad emission 
between 200 and about 600 kms" 1 . The velocity compo- 
nent at ~ 550 kms" 1 becomes weaker with rising rota- 
tional number. Emission of [C i] 1-0 is restricted to 200 and 
~ 500 kms" 1 only. In Circinus, the [Ci] peak temperature 
is a factor ~ 3 stronger than CO 4-3. 



Table 2. Observed line intensities in Circinus and NGC 4945. 
Calibration uncertainties are estimated to be 15%. The spec- 
tra in Figure[2]and the integrated intensities correspond to the 
listed angular resolution (FWHM). References: ICurran et al.1 
(1200131 ) .'' IMauersberger et all l|1996l ). 





Circinus 


NGC 4945 




line transition 


lint 


/int 


FWHM 




[Kkms" 1 ] 


[Kkms" 1 ] 


["] 


CO l-0 a 


180 


510 


45 


CO 2-1" 


177 


390 


38 


CO 3-2 6 




330 


38 


CO 3-2 6 


230 




15 


CO 4-3 


58 


212 


38 


13 CO l-0 a 


12 


30 


45 


13 CO 2-l a 


19 


45 


38 


[CI] 1-0 


163 


248 


38 



4. Physical conditions 

4.1. LTE 

In the optically thin limit, the integrated intensities of [Ci] 
and 13 CO listed in Tablc[2]arc proportional to the total col- 
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Fig. 2. CO 4-3 and [Ci] 1-0 spectra of Circinus and 
NGC 4945 obtained with the N ANTEN2 tele s cope togethe r 
with low- J CO obse r vation s bv lCurran et all (|2001aHl998h : 
iMauersberger et all (|l996t ) . All data are at 38" resolution, 
only the CO 3-2 spectrum of Circinus is at 15" resolution. 
The small amount of excess emission in the blue wing of 
the [C i] 1-0 line in Circinus maybe caused by small point- 
ing variations in the eastern direction as emiss ion in lower- J 
CO li nes extends to this region and velocities (jCurran et alj 
11998ft . 



4.2. Radiative transfer analysis of CO and 3 CO 

We modeled the 12 CO and 13 CO emission lines using an 
escape probability radiative transfer model for spherical 
clumps (IStutzki fc Winnewisserl I1985T) using the CO col- 
lision rates of ISchinke et al.1 (jl985l ). This non-LTE model 
assumes a uniform density and temperature in a homoge- 
nous clump. The physical parameters kinetic temperature 
Tkin, molecular density n(H2), and column density Aco de- 
termine the excitation conditions in this model (Tabk[3]). 

In NGC 4945 and Circinus, we used the ratios of the 
observed integrated intensities of CO 1-0 to 4-3 a nd the 
13 CO 1-0 and 2-1 transitions (|Curran et al.ll2001aft to ob- 
tain column densities, density and kinetic temperature, as- 
suming a constant 12 CO/ 13 CO abundance ratio. We used 
an a bundance rat i o of 40 in both sources in accordance 
with lCurran et all (|2001aft . The escape probility code uses 
an internal clump line width Av mo d which hardly effects the 
outcome of the model in the reasonable range of l-20kms _1 . 

In a simultaneous fit a x 2 -fitting routine then compared 
the J line ratios of the model output R J mod to the ratios of 
the observed integrated intensities R ohs and determined the 
model with the minimal \ 2 . We compute the normalized x 2 
with the degrees of freedom d = J —p and J being the num- 
ber of independent ratios and p the number of parameters, 
in our case Tki n , n(R2) and A*co, to be determined : 



mod obs 

3=1 



(1) 



umn densities. LTE column densities of carbon are rather 
independent of the assumed excitation temperatures (e.g. 
IFrerking et "aTlll989ft . We find N c = 3.4 - 3.9 10 18 cm- 2 in 
NGC 4945 and N c = 2.2 - 2.5 10 18 cnT 2 in Circinus for a 
temperature range of T ox = 20 — 150 K. 

We used the 13 CO J=l-0 and J=2-l integrated inten- 
sities to derive total CO column densities, assuming LTE, 
optically thin 13 CO emission, a CO/ 13 CO abundance ra- 
tio of 40 (jCurran et al.ll2001aft . and T GX = 20 K. We find 
a total CO column density of N C o = 1.0-1.7 10 18 cm~ 2 and 
N co = 4.1-6.7 10 17 cm- 2 in NGC 4945 and Circinus respec- 
tively depending on which transition is used, 1-0 or 2-1. 
Varying the temperature to T cx = 150 K the column den- 
sity of CO slightly increases up to Nqo= 3.2 10 18 cm~ 2 and 
N co = 1.410 18 cm" 2 for NGC4945 and Circinus. 

Another method uses the Galactic CO 1-0 to H2 
conversion factor Xmw = 2.3 1 20 cm ~ 2 (Kkms~ 1 )~ 1 
(jStrong et al.lll988l ; IStrong fc Mattoxlll996D and the canon - 
ical CO to H 2 abundance of 8.5 10" 5 (IFrerking et al.lll982f ) 
to derive Nco= 3.5 10 18 cm -2 for Circinus, i.e. a factor of 
~ 10 larger than the LTE estimate indicating that the X- 
factor is only 1/10 Galactic. As the abundance of CO maybe 
different in NGC4945 and Circinus this result has to be 
taken with caution. 

For NGC4945 IWang et all (|2004ft derive an X-factor 7 
times smaller than the Galactic value, which leads to Nco~ 
1.410 18 cm~ 2 , in good agreement with the LTE approxi- 
mation from 13 CO. The LTE column density derived from 
13 CO in Circinus also indicates an X-factor around 10 times 
smaller than the Galactic value. 

The CO/C abundance ratio is 0.29-0.50 in NGC 4945 
and 0.19-0.27 in Circinus using the LTE column densities. 



The errors crj due to calibration uncertainties are esti- 
mated to be 20%. 

In Circinus the 12 CO 3-2 line is missing, leaving one 
degree of freedom compared to two in NGC 4945. 

Tkin and n(H2) are determined with this step. To com- 
pare the modeled integrated intensities I mo d to the absolute 
observed intensities 7 b s we have to account for the veloc- 
ity filling, due to the velocity width Av mo( j of an individual 
clump to the width of the galaxy spectrum Av D b s and the 
beam dilution, due to the size of the modeled clump A c i 
compared to the beam area Abeam- 

The large velocity width of the observed spectra im- 
plies several clumps in the beam N c \=n Av b s /Av moc i with 
n > 1 (Table[3]). The beam dilution is determined by the 
fraction of modeled clump area to the beam size which we 
express in terms of an area filling factor per clump a ci= 
Ad/A beam- The total area filling factor is a — N c i 0a, ci- 
The size of the clump with radius R, A c i = irR 2 , can be 
inferred via the mass M and density ?i(H2) of the clump: 
R = (3/(47r)M/n) < - 1 / 3 '. In summary, the modeled inten- 
sities of the individual clumps / mo d arc converted to the 
intensities of a clump ensemble J Gns which can then be com- 
pared with the observed intensities: 



X 



(2) 



4.2.1. NGC 4945 

Figurc[5b shows the observed intensities of CO, 13 CO, and 
[C 1] together with two representative solutions of the ra- 
diative transfer calculations (see also Table[3|). 
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Table 3. Escape probability model results with two represen- 
tative solutions for each source. JVco and A r H 2 denote the total 
beam averaged column density of CO and H2, i.e. local clump 
column densities weighted by the total area filling factor <^a- 
M denotes the total mass. N c \ is the number of clumps in 
the beam, Aw mo( j is the modeled velocity width and </>a is 
the total filling factor to convert from modeled to observed 
intensities. The total cooling intensities are given in units 
of 10 _5 crgs~ 1 cm~ 2 sr _1 . 



Circinus 



NGC 4945 



x 2 


2.0 


9.6 


4.8 


12.4 


n(H 2 )ioc[cm 3 ] 


10 4 


10 3 


310 4 


10 3 


T kin [K] 


20 


100 


20 


100 


JVco [10 16 cm- 2 ] 


35 


50 


76 


63 


N c [10 16 cm- 2 ] 


230 


30 


330 


98 


jV H2 [10 20 cm- 2 ] 


37 


46.5 


89 


74 


M [10 6 M Q ] 


630 


792 


1385 


1114 


A^ m od [kms -1 ] 


5 


5 


10 


10 


Av b s [kms -1 ] 


186 


186 


188 


188 


N cl 


50 


38 


35 


40 


<t>A 


2.0 


6.3 


1.5 


4.0 


12 CO/ 13 CO abundance ratio 


40 


40 


40 


40 


CO/C abundance ratio 


0.15 


1.67 


0.23 


0.64 


[C 1] cooling intensity 


4.1 


7.88 


7.2 


11.8 


CO cooling intensity 


2.1 


2.8 


6.6 


5.7 


[C 1] / CO cooling intensity ratio 


2.0 


2.8 


1.1 


2.1 



CO. Fitting the CO and 13 CO lines, we find a degener- 
acy in the n(H2)-Tkin plane of the solutions for a rather 
constant pressure n(H 2 ) x Tki„ ~ 10 5 Kcm~ 3 . The best 
fits are achieved for a 12 CO/ 13 C O abundance ratio o f 
40, similar to the value found by ICurran et al.1 (|2001alh 
Low x 2 " va l ues ( sec Table© constrain the densities to 
"■(H2) = 10 3 — 10 5 cm -3 and temperatures to a wide range 
of Tki rl =20-180K with higher temperature solutions corre- 
sponding to lower densities. The best fitting solution is 
™(H2) = 3 10 4 cm -3 and Tki n =20K. However, this solution 
is not significantly better than e.g. n(H2) = 10 3 cm -3 and 
T kin =100K (Figib, Table©. 

The fitted column density Nco agrees with the LTE 
approximation for the CO-column density within a factor 
of 2-3. 

The peak of the modeled CO cooling curve at J=4 con- 
tains 35.6% of the total 12 CO cooling intensity of 6.6 10~ 5 
ergs _1 cm _2 sr _1 computed by summing the cooling in- 
tensity of the 12 CO transitions from J=l to 20 for the 
T kin =20K solution. 

We use the radiative transfer model to predict the 12 CO 
7-6 intensity. It is rather weak, depending strongly on T^m- 
It varies from k 1.3 Kkms -1 for the Tki n =20K solution to 
4.1 Kkms" 1 for t he Ikin=100K fit. 

ICurran et all (|2001aD find in this source 
n(H 2 )=3xl0 3 cm- 3 and T kin =100K from 12 CO obser- 
vations of the 3 lowest transitions and 13 CO data of the 



two lowest trans itions. In their multi-transition study 
IWang et all (|2004D estimate a density of n(H 2 )= 10 3 cm~ 3 
for an assumed temperature Tkin=50K from CO tran- 
sitions up to J = 3. In contrast, the observed CN and 
CH3OH lines indicate local densities around 10 4 cm -3 . The 
solutions for CO and 13 CO found in the literature are in 
good agreement with our parameter space of solutions 



showing that the additional CO 4-3 line does not help 
significantly to improve the fits. 

Atomic carbon. Assuming the same density, kinetic temper- 
ature, velocity filling, and beam dilution for carbon as for 
CO, we use the observed [Ci] intensity and the radiative 
transfer model to estimate the carbon column density and 
hence the CO/C abundance ratio. The CO/C abundance 
ratio varies between 0.23 and 0.64 for the two solutions 
listed in Table[3]The corresponding carbon column densi- 
ties arc N c = 3.3 10 18 cm- 2 and N c = 9.8 10 17 cm~ 2 for the 
high and the low density solution, rcspcctivly. The column 
density for the latter solution is a factor of 3 lower compared 
to the LTE carbon column density. The assumption of opti- 
cally thin emission in the LTE might be overestimating the 
column density compared to the escape probability mod- 
elling. We also use the model to predict the [C 1] 2-1 inten- 
sity. It is 16 Kkms" 1 for the T k in=20K/n(H 2 ) = 3 10 4 c m - 3 
solution. However, this result critically depends on the ki- 
netic temperature. The T kiri =100K/?i(H2) = 10 3 cm~ 3 so- 
lution yields 182Kkms- 1 . Also, the [Ci] 2-1/ [Ci] 1-0 line 
ratios change from 0.07 to 0.75, depending on the solu- 
tio.This shows that the high-lying transitions can be used 
to break the degeneracy. 

The total cooling intensity of both [C 1] lines is listed in 
Table[3]for both presented solutions. The cooling intensity 
of the two [C 1] lines is thus of the order of the total cool- 
ing intensity of CO with the C/CO cooling intensity ratio 
varying between 1.1-2.1 for the discussed solutions. 
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Fig. 3. Radiative transfer modeling results and observa- 
tions: Filled blue points show the observed CO and blue 
circles the modeled CO. The dotted line indicates the best 
fit solution. Filled red triangles show the observed 13 CO and 
red triangles the modeled 13 CO. Dash-dotted lines show a 
higher temperature solution. Filled black points are the ob- 
served [CI] integrated intensities and black circles the pre- 
dicted [CI] integrated intensities. 



4.2.2. Circinus 

CO. The modeled CO cooling curves are shown in Figure^ 
and the fit results are summarized in Tablc[3jThe CO line 
ratios are given for a A t; m od=5kms _1 so there are about 
40 clouds in the beam to account for the observed veloc- 
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ity width of Au o6s =186kms- 1 . The assumed 12 CO/ 13 CO 
abundance r atio of 40 is slightly lo wer than the values of ~ 
60, found bv lCurran etlll pOOlal) . 

Good fits corresponding to low x 2 (see Table© can be 
found for densities of n(H2) = 10 3 — 10 45 cm -3 and a large 
range of temperatures of T k i n =20-160K while the product of 
n(H.2) x Tkin stays approximately constant at ~ 10 s Kcm -3 . 
Again, a number of solutions provide consistent CO cooling 
curves. 

The lowest \ 2 is obtained for ra(H2) = 10 4 cm~ 3 , 
Tki n =20K and a column density of Aco = 3.5 10 17 cm -2 
assuming 50 modeled clumps in the beam. The column 
density Aco is well determined showing a steep gradient 
of x 2 -values for varying densities and temperatures. This is 
in reasonable agreement with the LTE-approximation from 
13 CO. A second solution at n = 10 3 cm~ 3 and T kin =100K 
also lies within the la-contour of the x 2 distribution (cf. 
Fig[3]and Tab©. 

The results agre e well with the solutions found by 
ICurran et ail (|2001a| ). They find T kin =50-80K and n(H 2 ) = 
2 10 3 cm -3 from observations of the 3 lowest 12 CO transi- 
tions and the 2 lowest 13 CO transitions. 

The modeled CO cooling curve peaks at J=4 contain- 
ing 35% of the total 12 CO cooling intensity of 2.1 10~ 5 
ergs- 1 cm- 2 sr" 1 for the T kin =20K fit. 

The predicted integrated intensity of 12 CO 7-6 is very 
weak, varying strongly from 0.1 Kkms -1 for the T k ; n =20K 
solution to 1.9 Kkms" 1 for T kin = 100 K. 

Atomic carbon. Assuming the same density and tempera- 
ture for carbon as for CO, we use the best fitting CO model 
to derive the [C i] 1-0 intensity, carbon column densities, 
and CO/C abundance ratios. The predicted CO/C abun- 
dance ratio is 0.15, again consistent with the optically thin 
LTE result. The predicted [Ci] 2-1 integrated intensity is 
6.1 Kkms -1 . 

As in NGC4945, changing the temperature has a large 
effect on the [Ci] 2-1 intensity. The T kin =100K solution 
yields 111 Kkms -1 and a higher CO/C abundance ratio 
of 1.67. Thus the [Ci]2-1/[Ci]1-0 line ratios changes from 
0.04 to 0.74. 

The corresponding carbon column densities for the two 
presented solutions are Ac = 2.3 10 18 cm -2 and Ac = 

3.0 10 17 cm -2 , repectivly. For the latter solution the column 
density is about a factor of 10 lower compared to the LTE 
carbon column density. The optically thin assumption for 
the LTE modelling is obviously not appropriate in the high 
temperature and low density scenario. 

The total cooling intensity ratio of [Ci]CO varies from 

2.1 to 2.8 for the presented solutions. Carbon is a stronger 
coolant than CO by a factor of 2-3. 

5. Discussion 

5.1. [Ci] 1-0 luminosities 

The [C i] 1-0 luminosities for the centers of the Seyfert 
galaxies NGC 4945 and Circinus are 91 and 67Kkms~ 1 kpc 2 
(Fig. ©. To date, about 30 galactic nuclei have been stud- 
ied in the 1- line of atomic c arbon , most of which are 
presen ted in iGerin fe Phillips! (|2000j ) and llsrael fc Baasl 
i|2002f) . The [Ci] luminosity of a source is an impor- 
tant property since it gives the amount of energy emit- 



ted per time which is proportional to the number of emit- 
ting atoms, i.e. proportional to the [Ci] column density 
in the limit of optically thin emission. The NANTEN2 
38" beam achieves ^700 pc resolution in Circinus and 
NGC 4945 which both lie at ~4Mpc. To achieve the 
same spatial resolution for galaxies at ^12Mpc dis- 
tance, e.g. for NGC 278 NGC 660, NGC 106 8, NGC 3079 
and NGC 7331 listed in llsrael fc Baasl (|2002fl . one would 
need an angular resolution of ~ 13", comparable to 
the 10" JCMT beam at 492 GHz. The luminosities stud- 
ied in these 7 sources thus all sample the innermost 
^700 pc. Area integrated [Ci] luminosities are found to 
vary strongly between ~l and ~ 160Kkms~ 1 kpc 2 i n 
these 7 galaxies (Fig.© (|Israelll2005t llsrael fc Baasl I2002D . 
Quiescent centers show modest luminosities 1 < L([Cl])< 
5Kkms _1 kpc 2 , while starburst nuclei in general show 
higher luminosities. The largest luminosities are found in 
the active nuclei of NGC 1068 and NGC 3079 which show 
50 and leOKkms^kpc 2 ([Israel fc Baasl I2002T ). NGC4945 
and Circinus also fall in this category (Fig.©. 

5.2. [Ci] 1-0 /CO 4-3 line ratios 

The [Ci] 1-0/CO 4-3 ratio of integrated intensities is 1.2 
in NGC 4945 and 2.8 in Circinus. For Circinus the ratio is 
larger than any ratios previously observed in other galac- 
tic nuclei or in the Milky Way. The [Ci] 1-0/CO 4-3 ratio 
is shown in Figure^] versu s the area integrated [C i] lumi- 
nosity. As also discussed in llsrael fc Baasl ([20021 ). we see no 
functional dependence. Galactic sources (not shown in this 
figure) would lie in the lower left corner, llsrael d2005h stud- 
ied 13 galactic nuclei and found that the [Ci] 1-0 line is 
in general weaker than the CO 4-3 line, but not by much. 
Ratios vary over one order of magnitude between 0.1 in 
Maffei2 and 1.2 NGC 4826. Galactic star forming regions 
like W3Main or the Carina clouds show much lower values, 
between about 0.1 and 0.5, which are consist ent with emis- 
sion from photon dominated regions (PDRs) dKramer et al. 
12004 I Jakob et al.ll2007t Kramer et aLll2007[ ). iFixsen et al. 
( 19991 ) find 0.22 in the Galactic center and 0.31 in the Inner 
Galaxy. The variation seen in the various Galactic and ex- 
tragalactic sources appears to be intrinsic and not due to 
observations of different angular resolutions. This is be- 
cause the frequencies of the two lines are very close and 
angular resolutions are therefore similar if the same tele- 
scope is used. 

5.3. [Ci] 1-0/ 13 CO 2-1 line ratios 

The [Ci] 1-0/ 13 CO 2-1 line ratios in NGC 4945 
and Circinus are 5.51 and 8.57, respectively (Fig.©). In 
NGC4945, the observed ratio is consistent wi t h re- 
sults found in previ ous studies (|Israel fc Baasl 120021 ; 
IGerin fc Phillipsl |2000D ranging up to ratios of 5. For 
Circinus, the ratio is again higher than any previous mea- 
surement. 

Two thirds of the sample of galaxies studied by 
llsrael fc Baasl ([20021 ) show [Ci] 1-0/ 13 CO 2-1 line ratios 
well above unity. The sample consists of quiescent, star 
burst and active nuclei. The highest [Ci] 1-0/ 13 CO 2- 
1 ratios are found in star burst and active nuclei, consis- 
tent with our observations. IGerin fc Phillipsl (|2000t ) find a 
similar result with two thirds of the galaxies in their sam- 
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Fig. 4. a [Ci]l-0/ 13 CO 2-1 ratios versus L([Cl]). b) [Ci]l- 
0/C O 4-3 ratios ver s us L([ C i]). Black points are the results 
from llsrael fc Baasl (|2002i ) and red points are the ratios 
from this paper. 



pie exceeding a ratio of 2. High ratios can be qualitatively 
understood in low column density enviroments with mild 
UV radiation fields. In these regimes most CO will be dis- 
sociated and the ga s-phase carbon will be neutral atomic 
([Israel fc B aas 2002?). This implies that dense, star forming 
molecular cloud cores are not the major emission source in 
galaxy centers. 

In general the studied centers of external galaxies 
show stronger [C i] emission than one would expect from 
Galactic observations, which show typical ratios of 0.2-1.1 
(jMookeriea et alj [2006) . In Galactic sources high ratios are 
found in low gas column densitis and medium UV radia- 
tion enviroments where 13 CO will be dissociated and atomic 
carbon remains neutral in the gas pha se i.e. in translucent 
clouds and at cloud edges (|Israelll20"05[ h He concludes that 
the dominant emission from galaxy centers does not stem 
fro m PDRs. 

iMeiierink et alj (|2007| ) studied irradiated dense gas 
in galaxy nuclei using a grid of XDR and PDR models. 
For the same density the predicted [Ci] 1-0/ 13 CO 2-1 
line ratios are significantly higher for the XDR- compared 
to PDR-models (Fig. 10 in Meijerink et al. 2007). We 
observed ratios of 61 and 95 in NGC4945 and Circinus 
respectivly, on the erg-scale. These ratios can be explained 
by XDR-models at high densities n(H 2 ) = 2 10 3 — 10 5 cm -3 . 
PDR-models explain the observed ratios in low density 
regimes with n(H 2 ) = 2 10 2 - 610 2 cm~ 3 . The high 
[Ci] 1-0/ 13 CO 2-1 line ratios observed in NGC4945 
and Circinus may hint at a significant role of X-ray 
heating in these galaxy nuclei as our predicted densities 
of n(H2) = 10 3 — 10 4 cm~ 3 agree w ith the high-density 
XDR-models in IMeiierink et all (|2007l ). 



5.4. Total CO and [Ci] cooling intensity 

In all sources studied bv lBavet et all (|200l but NGC 6946, 
the CO cooling intensity exceeds that of atomic carbon. 
They find the [C i] to CO cooling ratio to vary between 
~ 0.3 for M83 to - 2 for NGC 6946. NGC 4945 and Circinus 
show similarly high values as the latter galaxy: In NGC 4945 
this ratio is ~ 1 — 2 and in Circinus we find ~ 2 — 3. 



These values are around 2 orders of magnitude higher than 
the typical values found in Galactic s tar forming regions 
dJakob et alj|2007t iKramer et al.ll2007t) . 



5.5. Shape of the CO cooling curve 

The modeled CO cooling curve of NGC 4945 and Circinus 
peaks at J = 4. Observations of the higher lying CO lines 
i.e. the CO 6-5 and CO 7-6 lines will however be important 
to verify our model predictions and the importance of CO 
cooling relative to C. 

The shape and maximum of the cooling curve of 12 CO 
has been studied in a number of nearby and high-z galaxies. 
iFixsen et al.l (|1999( ) find a peak at J=5 in the central part 
of the Milky Way usin g FIRAS/COBE da ta and rotational 
transitions up to 8-7. iBavet et alj (|2006l ) observed 13 nu- 
clei in mid-J CO lines up to 7-6 and find that the peak 
of the cooling curves vary with nuclear activity. While nor- 
mal nuclei exhibit peaks near J up = 4 or 5, active nuclei 
show a peak near J up = 6 or 7. NGC 253 was observed 
at AP EX in CO upto J U p = 7 and shows a maximum at 
J = 6 (|Giisten et al.l [2006). On the other hand, studies of 
high rcdshift galaxies s how cooling curve s peaking as low 
as J=4 for SMM 16359 dWeifi et alJl2005D. as high as J=9 
in the case of the QSO APM 08279 (jWeiss et al.ll2007D . and 
peakin g at J=7 for the very high redshift (z = 6.4) QSO 
J1148 (IWalter et al1l2003h . 



5.6. Pressure of the molecular gas 

A wide range of temperatures Tk; n = 25 — 150K and den- 
sities n(H2) = 5 10 2 — 7 10 5 cm~ 3 has been found in sim- 
ilar studies of external nuclei, including ULIRGs, normal 
spirals, star burst galaxies, and i ntera cting galaxies (e.g. 
IBavet et all 120061: llsrael fc Baasl 120021) . In the irregular 
galaxy IC 10 rather high densities n(H2) ~ 10 6 cm -3 and 
low temperatures of 7km = 25K arc found while the molec- 
ular gas in the center of the spiral galaxy NGC 6946 is found 
to be l ess dense, nfHg ) = 10 3 cm~ 3 , but much hotter, Tkm = 
130K. iGiisten et ail (|2006l) studied NGC 253 to obtain 
n(H 2 ) = 10 3 - 9 cm" 3 and T kin = 60K while iBradford et~afl 
(|2003l ) investigated the same source and derived a higher 
T kin = 120K and density n(H 2 ) = 4.5 10 4 cur 3 . Both stud- 
ies used 12 CO 7-6 observations. However, the tempera- 
ture/density degeneracy cannot be resolved. The solutions 
we present in our study of NGC 4945 and Circinus show 
densities and temperatures of n(H2) = 10 3 — 10 4 cm~ 3 with 
a less well constrained temperature T^ n = 20 — 100K; de- 
pending on the density, as discussed before. 



5.7. [Cl] 2-1 /1-0 line ratio 

The modeled [Ci] 2-1/1-0 line ratios change from 0.07 to 
0.75 in NGC 4945 and from 0.04 to 0.74 in Circinus for 
the presented escape probability s olutions. Observed ra - 
tios vary from 0.48 in G333.0-0.4 (jTieftrunk et all |2001|) 
to 2.9 in W3 main for galacti c and extraga l actic sources 
(|Kramer et a l.l I2004D. For M82. IStutzki et all (|l998l ) found 
a ratio of 0.96. IBavet et al. I (I2006T observed ratios ranging 
between 1.2 (in NGC 253) to 3.2 (in IC342). The [Ci] 2- 
1/1-0 line ratio for the low temperature solution predicted 
for NGC4945 and Circinus is significantly lower than pre- 
vious results in the literature. 
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5.8. CO/C abundance ratio 

Compared to an a bundance ratio CO/C of 3-5 in NGC253 
(iBavet et al.ll2004D an d an average value of 2 in the nucleus 
of M83 ([White et alJ[l99l llsrael fc Baasll200ll ). NGC4945 
and Circinus have a higher fraction of atomic carbon in 
their nuclei which resembles the low CO/C abundance ra- 
tios found in Galactic translucent clouds. We find values 
of 0.23-0.64 in NGC4945 and 0.15-1.67 in Circinus. Due 
to the degeneracy of n(H2),T the abundance ratios can- 
not be determined more accurately. In galactic molecular 
clouds values range bet ween 0.16-100, low values are found 
in translucent clouds ([Stark fc van Dishoeckl Il994h while 
massive star forming re gions show high CO/C abundances 
(|Mookeriea et al.ll2006h . 

5.9. Future observations 

Future observations of CO 7-6 and [Ci] 2-1 will be im- 
portant to better understand the the kinetic temperature 
and density to resolve their degeneracy, and for an under- 
standing of the dominating heating mechanism i.e. X-ray 
or UV heating. We estimate detection of the [C i] 2-1 line 
with NANTEN2 within less than 10 minutes total observa- 
tion time under average weather conditions . The CO 7-6 
will be harder to detect, depending on the excitation con- 
ditions. At 810 GHz the 1 GHz bandwidth of the receiver 
translates into a velocity range of 370kms~ 1 . We therefore 
plan to stack future observations to be able to cover the 
broad line widths of Circinus and NGC4945. 
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